Is there an intermediate massive black hole in the Galactic center: 
Imprints on the stellar tidal-disruption rate 

Xian Chen\ k F. K. Liu^ 
ABSTRACT 

It has been suggested that an intermediate-massive black hole (IMBH) with 
mass 10^"^ Mq could fall into the galactic center (GC) and form a massive black 
hole binary (MBHB) with the central supermassive black hole, but current ob- 
servations are not sensitive to constrain all mass and distance ranges. Motivated 
by the recent discovery that MBHBs could enhance the rate of tidal-disruption 
events (TDEs) of stellar objects, we investigate the prospect of using stellar- 
disruption rate to probe IMBHs in the GC. We incorporated the perturbation 
by an IMBH into the loss-cone theory and calculated the stellar-disruption rates 
in the GC. We found that an IMBH heavier than 2000 Mq could distinguish- 
ably enhance the stellar-disruption rate. By comparing observations of Sgr A* 
with the fall-back model for stellar debris, we suggested that the TDE rate in 
our Galaxy should not significantly exceed 0.002 yr~^, therefore a fraction of the 
parameter space for the IMBH, concentrating at the high-mass end, can already 
be excluded. To derive constraint in the remaining parameter space, it is crucial 
to observationally confirm or reject the stellar-disruption rate between 10~^ and 
10"^ yr"^, and we discussed possible strategies to make such measurements. 

Subject headings: black hole physics - Galaxy: center - methods: analytical - 
X-rays: bursts 



Introduction 



It has been suggested that in the center of our own Galaxy, besides an unambiguously 
dete cted supermassive black hole (SMBH) of 4.3 x 10^ Mq whose location coincides with Sgr 
A* (IGenzel et al.ll2010l . and references therein), there might be another intermediate massive 
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black hole (IMBH), which is brought in by an infalling satellite galaxy (ILang et al.ll201ll ) or 
core-collapse star cluster ( iPortegies Zwart et al.ll2006[ ). The existence of massive black hole 
binary (MBHB) in the Galactic center (GC) would ease the tension between stellar dynamics 
and some observations, e.g. the formation of a young stellar disk inside 10" (about 4 pc), 
the central hole of size 5" discovered in the old stellar population, the apparent high velocity 
dispersion of a stellar association IRS 13E located at 3.5 " from Sgr A*, and the random 



orientation of the S-stars inside 1" (see Genzel et al. 2010. for a review). 



However, there has been no conclusive observational evidence for or against an IMBH 
at the galactic cen ter. Based on the p roximity of Sgr A* to the center of the Galac- 
tic nucl ear cluster (|Yu fc Tremaine and more exclusive ly on the proper motion of 



Sgr A* ( Hansen fc Milosavljevid 



2003 



iReid fc Brunthalerll2004l ). most parameter space can 
be excluded for IMBHs of masses Mimbh ^ 2 x 10'^ Mq residing at distances d > 2 
milliparsec (mpc) from Sgr A*. Slight improvement on the constraint can be obtained 
by taking into account the orbital stability of S-stars, which mostly excludes IMBHs of 
masses Mimbh ^ 10^ Mq at 0.3 mpc ^ d < 10 mpc (lYu fc Tremaind 120031). and 2 x 
10^ Mq < Mimbh < 4 x 10=^ Mq at 2 mpc < < 10 mpc JGualandris fc Merr"ittl I2OO9I ). 



Besides, the short coalescence timescale (< 10^ yr) due to gravitational wave radiation 
also argues against IMBHs o f any mass within a distance of ci ~ 0.1 mpc from Sgr A* 
(IHansen fc Milosavlievig 12003). Neverthe less, not all parameter space is excluded (e.g. see 
Figure 13 in iGualandris fc MerrittI boOQ ). The low-mass end at Mimbh < 2 x 10^ Mq is 
currently unaccessible, except that the IMBH should have d > 0.1 mpc to comply with the 
argument of coalescence timescale. Even in the high-mass end, it is not ruled out that an 
IMBH with 2 x 10^ Mq < Mimbh < 4 x 10^ Mq may reside at either 0.1 mpc <d<2 mpc 
or O.O6(Mimbh/1O^M0)2 pc < Mimbh < 3.2(Mimbh/1O^M0) pc. 

Previous studies on the stellar orbits bound to MBHBs discovered that the three-body 
interactions could significantly b oost the orbital ec centricities of the stars, partially due 
to the secular Lidov-Kozai effe ct (jivanov et al.ll2005l ) and more importantly due to chaotic 
interactions (jChen et al.ll2009l ). Consequently, compared to that in the single-black-hole 
case, a larger fraction of stars could reach the so called "tidal radius" of the central SMBH 



to be tidally disrupted, r esulting in a burst of stellar-disruption events (IChen et al.l 12011 
Wegg fc Nate Bodell201ll ). A stellar tidal-disruption event (TDK) will produce a powerfu l 



flare by releasing the gravitational energy of the stellar debris (jReesI Il988l : lUlmerl Il999l) 



and i s observable in mult i ple electromagnet ic bands ranging from radio to 7-rav (IKomossa 



2002 




Gezari et al. 


2009: 



Gezari et all l2009l : Ivan Velzen et all I2OIII : iBloom et all I2OIII : iBurrows et al.l I2OIID . 



Therefore, if IMBHs indeed reside in the GC, they should as well enhance the rate of TDEs, 
which may be imprinted in the emissions from the GC. 
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In this paper we first study the stellar-disruption rate in the GC in the presence of an 
IMBH. Based on our results, we discuss the prospect of constraining the parameters of the 
hypothetical IMBH using the observational TDE rates. The outline of the paper is as follows. 
In Section [2] we develop the loss-cone theory to incorporate the perturbation by an IMBH, 
then in Section [3] we calculate the stellar-disruption rate in the GC including the effects of 
the IMBH. In Section HI we investigate the observability of a TDE from GC based on the 
standard fall-back model for the stellar debris. Finally, we discuss our results in Section |5l 



2. Loss-cone theory 



Around the SMBH in our own Galaxy, the stars with pericenter distances smaller than 



rt 



5.4 



10"°-^ pc 



M. 



\ 1/3 







1/3 



(2) 



are subject to tidal disruption (lHillslll975t lReeslll988l ). where and m^, are, respectively, the 
radius and mass of star, and Rq and Mq refer to the solar values. For m^, = Mq and = Rq, 
this "tidal radius" is about 10 times greater than the Schwarzschild radius of the Galactic 
SMBH and has an angular size of 0.1 mas when viewed from the Earth. Tidal disruption 
creates a stellar-deficient region in the phase space of specific binding energy S and specific 
angular momentum J of stars. This region is conventionally referred to as t he "loss cone" 
because of its cone-like geometry when the system is spherically symmetric (IFrank fc Rees 
19761 : iLightman fc Shapiro! 119771 ). In the following, we restrict our calculations to the stars 
bound to the SMBH and far from the tidal radius, because they dominate the tidal- disrupt ion 
rate. This corresponds to an energy range of < £^ <^ GM,/rt, where a 75 kms~^ is the 
one-dimensional stellar velocity dispersion in the GC (IGenzel et al.ll2010l ). These stars are 
orbiting the SMBH on near-Keplarian orbits with semi-major axis a = GM,/2S, and the 
greatest semi-major axis we consider is Cmax = GM,/2a'^ ^ 1.6 pc. Given S, the maximum 
angular momentum is Jc = GM,/ (28)^^"^ and the boundary of the loss cone by definition is 
at Jie ~ {2GM,nf/\ 

The rate of tidal disruption is determined by the rate of stars diffusing into the loss 
cone, which is normally more efficient in the J direction of the £ — J phase space (J = 
|J|). In the simplest case where the system has a single stellar population with mass m^, 
suppose the successive mutual scattering between a star with energy £ and the background 
stars on average induces an angular- momentum change Jd{£) during one stellar orbital 
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period P{£) = 2nGM,/{2S)^/^ ~ 1.6 x 10^(0/0.1 pc)^/^ ^^en the "two-body" relaxation 
timescale can be calculated with T2b{S) ~ PJ^/ Jj^. The dependence of on the square of 
angular momentum reflects that the two-body scattering is a random process, so sums 
up incoherently. Given the number of stars n{S)dS in the energy range S ^ S + dS, the 
loss-cone filling rate is proportional to n{£)d£ /T2h{£)- More careful analysis, taking into 
account the detailed distribution function f{S, J) at the loss-cone boundary and the fact 
that stars can be deflected into and out of the loss cone during one P{S) when Jd ^ Jy^, 
gives the following form for the loss-cone filling rate due to two-body relaxation: 



J'2bi£)d£ = 3ln{£)dE/P{£), 



(3) 



where 



(4) 



Perets et al.l 120071 ). To calculate {Jd/JicY, we adopted Equation (14d) in 



(Yound 


1977; 


Perets et al. 


2007) 



Bahcall fc Wolj f ll977h . 



For stars within about 0.1 pc from the Galactic SMBH, another efficient relaxation pro- 
cess is resonant relaxation (RR), a coherent chan ge of J driven by the torque exerted by the 
grain y gravitational field of finite number of stars ( iRauch fc Tremainelll996l : iRauch fc Ingalls 
19981). Our scheme to calcula te the loss-cone filling rate due to RR is analogous to that in 
Hopman &: Alexanderl ( l2006al ). but we also revised their formula to include the dep endence 
of RR on stellar orbital eccentricity e = [1 — (J/Jc)^]^/^ f Giirkan fc HopmanI 120071 ). In the 
standard case with a single black hole (BH), the timescale for coherent variation of J by RR 
is limited by the orbital precession timescales induced by general relativity (GR) and by the 
non-Keplerian potential of the surrounding stellar cusp. Given P{£) and e = (1 — e^)^/^ for 
a stellar orbit, we calculate tcR = (2e^/3) (a/r5)P(i^) for the GR precession timescale and 
tM — e~^{M,/m^)P{£)/N{> £) for the cusp-induced precession, where rs = 1GM,j(? is 
the Schwarzschild radius and A^(> £^ denotes the number of stars with energy greater than 
£. Then the coherent variation timescale, t^, can be derived from 1/t^ = |l/tGR — 1Am|, 
the minus sign before 1/tA/ due to the opposite precession directions. During t^, J varies 
coherently by AJ^j ^ JRRt^, where >/rr — 0.2 5eA^^/^('> £)Gm^:/a is the RR torque exerted 
on the stellar orbit ( iGiirkan fc Hopmanll2007l ). On longer timescale t ^ t^, the coherence is 
broken due to orbital precession, therefore J varies incoherently as A J ^ /S.J^it/t^Y^'^ . The 



RR timescale to erase the initial angular momentum, J 
and more precisely 



eJc, is then Trr ~ {€Jc/AJ^)\ 



Trr(^, J) 



2.55e2 P\£) fM, 

e2 iV(> £) 



(5) 
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( IGiirkaii fc Hopmanll2007l ). The average relaxation timescale Trr for the stars with the same 
energy E is given by the integration of the equation d.P / = dt/T^^{£, J), i.e., 



(6) 



Trr(^) = / TnK{S,J)djyjl 

Then loss-cone filling rate due to RR can be calculated with 

FKK{S)d8 = n{S)dS/fnn{S). 



(7) 



If an IMBH with mass Mimbh (the mass ratio of the IMBH-SMBH binary being q = 
Mimbh/M,) resides at a distance d from Sgr A*, the tidal force of the IMBH also exerts 
a torque on a stellar orbit, which induces an additional coherent variation of J at the rate 
Jk - Jc/Tk, where 



K 



— (-) 

3nq \dJ 

I6V2 /'a\l/2 



P(a) 



(a < d/2) 
(a > d/2) 



is analogous to the Lidov-Kozai timescale f Lidovlll962l : lKozailll962l ) but also acco unts for the 
stars in the chaotic regime with a ~ (i which dominates the loss-cone refilling (jChen et al. 



2OO9I . I2OIII ). We note that during chaotic interactions the stellar orbits could reach extreme 
eccentriciti es (e ^ 1) irrespe ctive of their initial inclinations or the 2:-components of angular 
momenta (jChen et al.ll201ll ). which are fundamentally different from those interactions in 
the secular Lidov-Kozai mechanism and significantly increase the stellar reservoir for tidal 
disruption. Because of the extra nodal precession induced on the stellar orbit by the IMBH, 
the coherent vari ation of J is limite d by the new timescale t'^ = \l/tK + 1/^gr ^ ^/t 



M\ 



where — eT^- (jivanov et al.ll2005[ ). During t'^, the variation of J^, due to both RR and 
the IMBH perturber, is A{J'^y = ( + The resulting timescale for erasing the 

initial J is Tco{S, J) ~ t^(eJc)^/A( J^)^. Therefore, the coherent relaxation timescale due to 
both RR and an IMBH is 



fun 



j2 



n,i£,j)djyj, 



c ' 



and the corresponding loss-cone filling rate is 



(9) 



(10) 



We assume /cj — 0.5 to correct, to the zeroth order, the star loss due to slings hot ejection, 



since /ej is a complex function of q, a/d, and the orbital eccentricity of the IMBH (jChen et al. 

201 ih . 
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3. Stellar-disruption rate in the GC 



T o give n{£)d£ and A^(> £), we adopted a stellar distribution for the GC from lSchodel et al. 
(120071 ). whose mass volume density has an outer Bachall-Wolf and an inner 7 = 1.2 power-law 
profile with a break radius at ~ 0.22 pc. The corresponding stellar distribution function 
f{S) scales as f{8) oc ^"P, where p = 1/4 (-0.3) for £ < GM,/n {£ > GM,/n). The 

'2{GM,/r - S)f{S)dS, then 



jGM./2rt Using these 



normalization for f{£) is derived from p(r) = 47rm* J^^^'^^ 
by definition we derive n{S) = 4:7r'^J^f{S)P{S) and N{S) 
quantities, as well as Equations ( fTOl) and (j3]), we calculated the stellar-disruption rates in 
the GC for a grid of hypothetical IMBHs of different Mimbh and d. 

The dashed line in Figure [1] shows the stellar-disruption rate contributed by two-body 
relaxation at different semi-major axis, assu ming m^, = 1 Mp^. Since Jp oc ml^'^ and r* oc m°'^ 



for main-sequence stars when m* < 20 Mq ( iKippenhahn fc Weigertlll990l ). the rate J^2h scales 
only mildly with m*, as 1/ Inm* in the limit Jd ^ ^ic and m^'^^ in the limit Jd ^ ^ic (see 
Equation (jlj)). We do not show the contribution from the stars at a > 1.6 pc, because 
the Keplerian assumption would break down for these stars, but their contribution to the 
integrated loss-cone filling rate N = J J^{S)dS is insignificant anyway. Inside the central 
0.004 pc, the stell ar mass is dominated by comp act objects such as white dwarfs and stellar- 
mass black holes (IHopman fc Alexander! l2006bl ). therefore the rates inside 0.004 pc are not 
shown either. The total stellar-disruption rate due to incoherent relaxation, integrated over 
the range 0.004 pc < a < 1.6 p c is A^^h ^ 3.5 x 10~^ yr"^, agrees well with the previous more 



sophisticated calculations (e.g. iMerritt Il2010l ). 



The dotted line in Figured] shows the contribution to stellar-disruption rate by RR only, 
without the perturbation of an IMBH. Compared t o two-body relaxation (dashed line), RR 
becomes more important at a < 0.1 pc (also see iHopman fc Alexanderl l2006a( ). and the 
sharp decline of J^{S)S inside a ~ 0.01 pc is caused by the quenching of RR by relativistic 
precession. The integrated stellar-disruption rate is A^rr 
depend on because both n{£)d£ and T'rr(£^) scales as 



3.5 X 10 ® yr ^, and it does not 



If an IMBH with Mimbh = 10^ Mq resides in the GC, the total stellar-disruption rates 
contributed by the coherent relaxations, i.e. RR plus IMBH perturbation, are shown in 
Figure [T] as the solid and dot-dashed lines. The rates are derived under the assumption 
= Mq and scale as if varies. Comparing the solid lines with the dotted one, we 
found that the presence of the IMBH dramatically enhances the loss-cone filling rate. The 
integrated stellar-disruption rate derived form the solid lines, A^co, increases with decreasing 
d when > 0.1 pc. It reaches a maximum of Aco — 1.1 x 10"'^ yr~^ when d ~ 0.07 pc, which 
is a factor of 300 greater than the unperturbed value of iVRR, and 30 times higher than A'2b. 
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10-1 




10"^ 10"i 10' 

a / pc 



Fig. 1. — Contribution to loss-cone filling rate by stars at different semi-major axis, assuming 
= Mq. Different stellar relaxation mechanisms are indicated with different line styles, 
i.e., dashed for two-body relaxation, dotted for RR only, and solid for RR+IMBH with 
^iMBH = 10^ Mq. The annotations on the lines are the assumed distances {d) in unit of 
pc for the IMBH. An IMBH with 10^ Mq at c/ = (2 x lO'^, IQ-^, 8 x IQ-^) pc is excluded 
by observations and thus the loss-cone filling rates due to the mechanisms RR+IMBH are 
shown in dot-dashed lines. 
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As d further decreases from d = 0.07 pc, the integrated rate A'^co becomes smaller, because 
the coherent variation of J at a ~ d becomes more susceptible to quench by relativistic 
precession. 

Figure |2] (intensity map and contours) shows the total mass disruption rate, M = 
'm^.{N2h + Aco), which is insensitive to the assumption of m^,, as a function of Mimbh and d. 
The dashed line indica tes the parameter space tha t is excluded by the observed dynamics 
of S-stars and Sgr A* ( iGualandris fc MerrittI 120091 ). In general, M is a increasing function 
of MjMBH (because N^o oc Mimbh), and given Mimbh, the rate peaks at ~ 0.01 — 0.5 pc. 
For example, when Mimbh = (10^ 10"^, 10^) M©, the maxima of M are (8.3, 110, 1600) x 
10~^ Mq yr~^, occurring at d = (0.03,0.07,0.3) pc. To the left-hand-side of the dashed 
line, where the possibility of an IMBH cannot be excluded by the current observations, the 
perturbed stellar-disruption rate ranges from 3 x 10~^ Mq yr~^ when Mimbh ^ 10^ Mq to 
as high as 4 X 10~^ Mq yr"^ when Mimbh — 10^ Mq. We also derived Nco for q = 1/81, 
1/243, and 1/729 from t he data of our prev ious scattering experiments in pseudo-Newtonian 
gravitational potentials ( IChen et al.ll201ll ). The resulting numerical rates agree with those 
from the above analytical calculations within a factor of 2, despite many simplifications in 
our analytical model. 



4. Observability of TDEs in the GC 

The black contour in Figure |2] (M ~ 10~^ Mq yr~^) indicates that an IMBH with 
Mimbh > 2000 Mq, if currently resides in the GC, would significa ntly raise the possibility 



of TDEs in the past 10^-10^ years. The standard model for TDK flReeslll988l : lUlmerl 119991 ) 
predicts that a fraction of the stellar debris will "fall back" to the periastron, forming an 
accretion disk and producing a shock- heated hot spot, whose initial luminosity should be 
close to or even exceeding the Eddington luminosity of the central SMBH [Le — 5.5 x 
]^q44 gj,g g-i GC). Is the hypothesized high stellar-disruption rate consistent with the 
observed quiescent state of Sgr A*? 

In the framework of the "fall back model" , the most bound seller debris has a binding 
energy of AS ~ kGM,r^/ {rt/ and will return to the periastron after a time delay of 

trr^in = 27rGM. (2 A^) "^/^ ^11) 

- Q.22k-^/^(5-\m,/MQ)-\r,/ R^f/^ yi (12) 



( iLi et al.ll2002[ ). where A; ~ 1 3 accounts for the stellar spin at periastron before tidal 

disruption and /3 = rj/rp ~ 1 denotes the ratio between the tidal radius and the pericenter 
distance Tp. About a fraction / ~ 0.5 of the stellar debris remains bound to the central SMBH 
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Fig. 2. — Integrated mass disruption rates (intensity map and contours), including contribu- 
tions from both coherent and incoherent relaxations, for a grid of Mtm rh and d. The param- 



eter s pace for IMBH excluded by the dynamics of S-stars and Sgr A* ( iGualandris fc Merritt 



20091 ) is the region to the right hand side of the dashed line. 
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( lEvans fc Kochaneklll989l ). and the resulting fall-back rate as a function of time is 



fb 



2fm^ ft — to 



X 



S^min 



-5/3 



1/3 



Rq 



0.12 Mq yr"^ fc"^/?"^ 
1 yr 



(13) 



where Id denotes the time of stellar disruption and t — to > tmm- We note that the fall- 
back rate when Mimbh > 4 x 10^ Mp, and d > 0.0 03 pc is not calculated self-consistently, 



because in this case, according to iLiu et al.l (120091 ) . Mfb would be interrupted at t — tn > 
10 yr. Anyway, this pararn eter space is excluded by the dynamics of Sgr A* and S-stars 
( IGualandris fc Merrittll2009[ ). Because of mutual collisions when returning to the periastron. 



the fall-back material circularizes at a radius of about 2rp, and the prompt release of the 
kinetic energy gives rise to a hot spot with bolometric luminosity 



Lhs = eMfbC^ ^ 9.5 x 10^^ erg s~^ 

2/3 



0.5 



Mrr 







t-tp 

1 yr 



-5/3 



(14) 



where e ~ GM,/ (4rpC^) is the conversion efficiency of kinetic energy into radiation and c is 
the speed of light. The spectral energy distribution (SED) of the hot spot is mostly likely 
black-body (BB) because of the large opacity of the fall-back material. The effective tem- 
perature of the BB depends on the emission area, Abb? which is uncertain due to the current 
poor understanding of the structures of the colliding streams flKochanekl Il994l : iKim et al. 
19991 ). Since the vertical scale hight orthogonal to the orbital plane of the stellar stream 
is greater than r*, the cross section for stream collision should not be much smaller than 
27rr^,rp, and the geometry factor C, = ABB/(4vrrp) for the emission area is likely greater than 
r*/rp ~ 0.01. For a typical value of C = 0.1, the effective temperature of the hot spot is 



T, 



BB 



L 



bol 



\asA 



BB 



1/4 



bol 



~ 1.8 X 10° K — 



(3 



k 



1/4 



/ 



0.5 



1/4 



1/4 / N 1/3 



0.1 



-1/4 / ^ \ -1/2 . , , 



1 yr 



Mo 

-5/12 



(15) 



where as is the Stefan-Boltzmann constant. 



To get a sense of the brightness of the hot spot in the GC, we calculated Mfb and Lhs as 



functions of t — tp with the fiducial parameters (/, k, /3, (^,r* 



(0.5,1 



,1,0.1, R, 



Me). 
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We also increased k and /3 to /c = 3 and /3 = 5 to estimate the lower limits to Mfb and 
Lhs- The results are given in Table [1] with the lower limits bracketed by parentheses. We 
found that 10 years after a TDE in the GC, the fall-back rate has already dropped to 
about 30-2600 times smaller than the Eddington rate, defined as Me = Le/{0.1c^) ~ 8.4 x 
10"^ Mq yr~^. By this time, the bolometric luminosity of the hot spot is 4 — 5 orders 
of magnitude smaller than the Eddington luminosity, when taking into account that e is 
much smaller than the typical value of 0.1 derived from the disk accretion theory. We note 
that when t — tu > 100 yr, the bolometric luminosity contributed by the accretion disk 
produced by the circularized material is ne gligible because the fall-back rate is too low t o 
support a radiatively efficient accretion disk (jCannizzo et al.lll990l : iMenou fc Quataertll200ll ) . 
Even in the e xtreme case that the accretion flow can instantly drain the fall-back material 
Juimei]ll999L ;he luminosity of the radiatively inefficient accretion flow" is still as low as 
IQ38 s-i(Mfb/10"^ Mq YT~^f jNaravan fc Yll995l : [Naravan et al1ll998l). Observations o f 
Sgr A* revealed that the bolometric luminosity is of order io^^~^^ erg s~^ (IGenzel et al.ll2010l ). 



and the SED c an be well modeled w ith a radiatively inefficient disk with an accretion rate of 
10"^ Mq yr~^ (INarayan et al.lll998l ). Detailed modeling of the accretion flow also indicates 
that the accretion rate inside 20 gravitational radii (smaller than the typ ical tidal radius) of 
the c entral SMBH is even smaller, ranging from 10"^ to 10"^ Mq yr~^ (jShcherbakov et al. 



2OI2I . and references therein). These facts, when compared to Mfb and Lhs derived from the 
fall-back model, suggest that TDEs are unlikely to have occurred in our own Galaxy within 
the last 4 — 5 centuries. 

We also calculated Tbb of the hot spot according to Equation (HM . The results for 
the flducial parameters (data without parentheses) and for {k, /3, () = (3, 5, 0.01) are also 
shown in Table [H According to Wein's displacement law, the SED of the hot spot initially 
peaks at the UV band and will shift toward IR as t — t^^ increases to 10^ years. Since 
IR bands have relatively low extinctions toward the GC and are commonly used by the 
ground-based telescopes to monitor the GC, we calculated the monochromatic flux of the 
hot spot at a wavelength of 2.2yum (i^'-band) with an extinction of A = 3. The results 
are given in the last two rows of Table [H where the data in parentheses correspond to 
non-flducial parameters with {k,P,Q = (3,5,0.01). For other non-flducial combinations of 
(A;, /3, C, /, r,,), one can use the scaling oc Ly^oi/T^s to derive the IR flux. The K- 
band flux from the radiatively inefficient accretion flow is very difficult to derive because of 
many unknown mode l parameters, but is typically smaller than 10(Mfb/10~^ Mq yr~^)^ mJy 
(INarayan et al.lll998l ). Observationally, the K- hand flux of Sgr A * is measured to be < 5 
mJy with occasional flares as bright as 30 mJy (IMorris et al.l 120121 . and references therein). 
Despite the uncertainties in the model predictions, we found the null hypothesis that a TDE 
occurred in the GC within the last 7 — 8 centuries to be in severe conflict with the observed 
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i^-band flux of Sgr A*. 



5. Discussion 

In order to test the possibility and constrain the parameters of a hypothetical IMBH 
in the GC, we studied the effects of an IMBH on the stellar-disruption rate in the frame- 
work of the loss-cone theory. We found that an IMBH with mass greater than 2000 Mq 
could significantly enhance the stellar disruption rate in GC (Figure [1]) because it coherently 
perturbs the angular momenta of stars in a fashion analogous to but more chaotic than 
the Lidov-Kozai effect. The maximum mass disruption rate increases with black hole mass 
as 10"'^ M4 Mq yr~^, where M4 = Mimbh/10^ Mq, and is reached when the IMBH arrives 
at about 0.1" ~ 1" from Sgr A* (Figure [2]). These results have already brought up some 
intriguing implications. For example, it is speculated that the high velocity dispersion of 
IRS 13E, a young stellar association locating at 0.14 p c from Sgr A*, may be induced by 



an embedding IMBH w ith mass greater than 10^ Mq flMaillard et all l2004t iSchodel et all 



2OO5I : iFritz et al. II2OIOI ). If so, the mass disruption rate in the GC according to our calcu- 
lation should be greater than about 5 x 10~^ Mq yr~^, up to about 0.02 Mq yr~^ in the 
extreme case with Mimbh = 10^ Mq. The results also imply that an IMBH in the GC would 
greatly increase the c hance for us to detect infalling objects on the course toward Sgr A* 



fiGillessen et al.ll2012h . 



Is the hypothesized high stellar-disruption rate consistent with the current quiescent 
state of Sgr A*? To estimate the set-off time for the most recent TDK in the GC, we 
calculated the evolution of the bolometric luminosity and the IR flux of a TDE in the GC 
using the fall-back model (Section H]). When comparing the model predictions with the 
observations of Sgr A*, we found that the current quiescent state of Sgr A* only allows 
a stellar-disruption rate lower than about 0.002 yr~^. According to this limit, the IMBH 
is unlikely to reside in the right-hand side of the log(M) = —2.7 contour in Figure |2l 
Interestingly, most part of this disfavored parameter space coincides with the space already 
excluded by the dynamics of Sgr A* and the S-stars (dashed line in Figure [2]), confirming 
the previous constraints on the IMBH. In addition, the limit on M excludes a patch of the 
parameter space at (i < 0.002 pc and M > 0.002 Mq yr~^, which was not reached by the 
dynamics of Sgr A* and the S-stars. Furthermore, to comply with the limit on M, the 
hypothetical IMBH in IRS 13E should not exceed 2 x 10^ Mq. 

The above success encourages us to investigate the prospects of using stellar-disruption 
rates to probe even smaller and closer (to Sgr A*) IMBHs. For the parameter space of IMBH 
that is still allowed by the dynamics of Sgr A* and the S-stars (Figure [2]), the maximum 
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stellar-disruption rate is M ~ 10"^ Mq yr"^ when Mimbh = 10^ Mq and M ~ 10""^ Mq yr"^ 
when Mimbh = 10^ Mq, and these maxima occur at ~ 0.04 pc (about 1 arcsec). Besides, 
at a distance as close as ci < 0.002 pc (about 50 mas), where there is nearly no constraint on 
Mimbh by the previous observations, the maximum M varies from about 4 x 10~^ Mq yr~^ for 
Mimbh = 10^ Mq to about 10"^ Mq yr'^ for Mimbh = 2000 Mq. Combining these results 
and taking into account that the stellar-disruption rate for a single SMBH in the GC is 
7V2b — 3.5 X 10^^ yr"^, we suggest that a stellar-disruption rate between 10~^ and 10~^ yr"^, 
if identified in our own Galaxy, would be a strong evidence for an IMBH more massive than 
2000 Mq lying in the central parsec of the GC. On the other hand, if the stellar-disruption 
rate higher than 10~^ yr~^ is rejected by future observations, a large parameter space, which 
is not covered by the current constraints on the hypothetical IMBH, could be excluded. 

A stellar-disruption rate of 10^^ yr^^ is already too low to be measured directly from 
decadal observational cycles, therefore measuring a rate down to 10~^ yr~^ requires special 
strategies. It is o bserved that the cen tral tens of parsecs of the GC are populated by molecular 
clouds (MCs; see lGenzel et al.ll2010l . for a review). A TDE at Sgr A* could leave imprints on 
the MC system in the form o f reflected hard X-rays or fluorescent iron lines decades or even 
centuries after the outburst (lYu et al.l l201ll ). A systematic search for such "light echoes" 
within 30 pc from Sgr A* indeed revealed a ma jor flare as lumino us as 3 x 10^^ erg s~^ in 
X-ray starting to decay at about 150 years ago (IMorris et al.ll2012l . and references therein). 
These observations give a upper limit of 0.004 ~ 0.006 yr~^ to the TDE rate in the GC 
because according to Equation f[T^ it takes at least 180 — 280 years (depending on k and (3) 
for the luminosity of hot spot to drop below 3 x 10'^^ erg s~^, not to mention that accounting 
extra X-ray contribution from an accretion disk would extend the time elapse t — to- To test 
whether the flare recurs every 200~300 years, it is crucial to look for a second light echo in the 
MC system 30 — 100 pc away from Sgr A*. On longer timescales ranging from thousands to 
millions of years, the imprints of TDEs could be detected in the forms of enhanced electron- 



posit r on an nihilation lines, high-energy gamma rays, and cosmic rays (jCheng et al.l 12006 



20071 . 12OIII ). If absorbed by the bio-circle, the records of en hanced high-energy emissions 



could be preserved on the earth surface for thousands of years flMiyake et al.ll2012l ). In order 
to constrain M in the GC over an extended timescale of 10^~^ years, we appeal for more 
efforts in measuring and modeling of the emissions in these high-energy bands. 
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Table 1. Predictions from the Fall-back Model 



t-to (years) 


10 


102 


10^ 


10^ 


Mfb (Mo yr-i) 


10-2.6 


10-4.3 


10-5.9 


lo-^-"^ 




(10-^-5 


10-6.1 


10-7.8 


10-9.3) 


Lhs (erg s-i) 


1041.3 


1039.6 


1038.0 


1036.3 




(1040.1 


1038.5 


1036.8 


1035.1) 


Tbb (lO^K) 


79 


31 


12 


4.5 




(160 


61 


23 


8.9) 


F,(i^) (mJy) 


160 


56 


15 


2.0 




(35 


12 


4.1 


0.96) 



